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Abstract 

Based on interaction potentials between a heavy quark and antiquark as extracted from recent QCD 
lattice calculations, we set up a Brueckner-type many-body scheme to study the properties of light (anti-) 
quarks in a Quark-Gluon Plasma at moderate temperatures, T ~ 1-2 T c . The quark- antiquark T-matrix, 
including both color-singlet and -octet channels, and corresponding quark self-energies and spectral functions 
are calculated self-consistently. The repulsive octet potential induces quasiparticle masses of up to 150 MeV, 
whereas the attractive color-singlet part generates resonance structures in the q-q T-matrix, which in turn 
lead to quasiparticle widths of ~200 MeV. This corresponds to scattering rates of ~1 fm _1 and may reflect 
liquid-like properties of the system. 
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I. INTRODUCTION 




A central goal of the relativistic heavy-ion collision program is the creation and identification 
of new forms of highly excited nuclear matter, in particular a deconfined and chirally symmetric 
Quark-Gluon Plasma (QGP). At sufficiently high temperature T, due to asymptotic freedom of 
Quantum Chromodynamics (QCD), the QGP is expected to be a weakly interacting gas of quark- 
and gluon-quasiparticles with comparatively small thermal masses, rriq t g ~ gT. Recent data from 
the Relativistic Heavy-Ion Collider (RHIC) indicate, however, that the produced matter exhibits 
strong collective behavior which is incompatible with a weakly interacting QGP: standard (2 «-> 2) 
perturbative QCD (pQCD) cross sections for quarks and gluons do not allow for rapid thermaliza- 
Jl| as required in hydrodynamic models to reproduce the observed magnitude of the elliptic 
, 0, 0] • With estimated initial energy densities well in excess of the critical one predicted by 
lattice QCD (1QCD), e c ~ 1 GeV/fm 3 , the question arises what the nature of the produced medium 
at temperatures T ~ 1-2 T c is (T c ~ 170 MeV: critical temperature). Of particular importance 
is the identification of the relevant interactions that can lead to sufficiently large scattering rates 
while maintaining consistency with the QGP equation of state (EoS), as determined in 1QCD. 

Recent (quenched) 1QCD calculations found intriguing evidence that mesonic correlation func- 
tions, after transformation into Minkowski space, exhibit resonance (or bound-state) like structures 
for temperatures up to ~2T C . This was first observed for low-lying charmonia (rj c , J/ip) 0, @, 0]> 
but subsequently also for mesonic systems with lighter quarks p, As is well known, resonance 
scattering is typically characterized by isotropic angular distributions and thus more efficient in 
randomizing momentum distributions than forward-dominated pQCD cross sections. Indeed, a re- 
cent calculation based on the assumption of resonant "D" -meson states in the QGP has shown 
that thermal relaxation times for charm quarks are reduced by a factor of ~3 as compared to using 
perturbative rescattering cross sections. The possibility of light hadronic states (especially for the 
pion and its chiral partner a) surviving above the phase transition has been suggested some time 
ago using effective quark interactions, e.g., within the Nambu-Jona-Lasinio model [ill. I12J . within 
the instanton-liquid model based on euclidean correlators jl3| , or more recently in Refs. jlJ. fl5| . 

To make closer contact to 1QCD, some recent works have extracted a (color-singlet) heavy-quark 
(Q-Q) potential, V\, from the corresponding 1QCD free ener gy, F\ , at finite T, and injected it into 
a Schrodinger equation to infer quarkonium properties Il7l . Il8j | . Reasonable consistency was 
found in that the heavy-quark bound states dissolve at roughly the same temperatures at which the 
peaks in the 1QCD spectral functions disappear (~2T C for J/i/j and rj c ), provided the free energy 
was converted into a potential by subtracting an entropy term according to V\ = F\ — T dF\jdT . 
A similar approach has also been applied to the light-quark sector in Refs. 

HI 0, EJ, where 

the q-q potentials from unquenched 1QCD (including colored channels) have been supplemented 
by relativistic (and instanton-induced) interaction corrections. Assuming rather large quark- and 
gluon-quasiparticle masses, m q g ~ 3-4 T (motivated by 1QCD calculations of temporal masses 22]), 
light mesonic, as well as a large number of colored diquark, quark-gluon and gluon-gluon, bound 
states have been found. Both quark-/gluon-quasiparticles and binary bound states together were 
shown to approximately reproduce the EoS from 1QCD. However, the effects of finite widths for 
both (anti-) quarks and bound states, which are essential to address scattering problems, were not 
included. 

In the present article we emplo y q uark-antiquark potentials extracted from 1QCD (including 
relativistic corrections as in Refs. |l9l . Eol . EH ]) within a 3-dimensionally reduced Bethe-Salpeter 
equation to evaluate (anti-) quark interactions in the QGP. We compute the pertinent scattering (T- 
) matrices in both color-singlet and -octet channels and calculate the quark self-energies including 
both real and imaginary parts (corresponding to quasiparticle masses and widths). The self- 
energies, in turn, are reinserted into the q-q propagator of the T-matrix equation, constituting a self- 
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consistency problem which we solve by numerical iteration. We comment on possible consequences 
of our results for quasiparticle masses and widths with respect to the QGP EoS and (anti-) quark 
rescattering timescales, respectively. 

Our article is organized as follows. In Sec. [Q] we present our parametrization of 1QCD data 
for the singlet free energy and extract a pertinent quark-antiquark potential including both color- 
singlet and -octet contributions. In Sec. Mil we set up our self-consistency problem comprising the 
q-q scattering equation and in-medium single particle self-energies and propagators, and discuss the 
underlying assumptions and approximations. The numerical results with accompanying discussion 
for the T-matrix and self-energy in a nonperturbative QGP are contained in Sec. IIVI In Sec. IVlwe 
conclude and give an outlook. 



II. QUARK-ANTIQUARK POTENTIAL FROM LATTICE QCD 

To obtain a driving term (potential) for a q-q scattering equation we take recourse to 1QCD 
calculations of the static free energy for a Q-Q pair. The Bielefeld group has performed extensive 
studies of this quantity based on Polyakov loop correlators for both the pure-glue SU(3) 24], 2f\ 
and iVy=2-QC p l2q.l27ll. Various parameterizations thereof have been given in the literature, cf., 
e.g., Refs. Q, EJHH . For the temperature range T = 1.1-2 T c , it turns out that unquenched 




FIG. 1: Left panel: lattice QCD results for the color-singlet free energy from unquenched simulations [2jj 
for 6 different values of the temperature (symbols) compared to our fit function, Eq. represented by the 
various curves. Right panel: corresponding potential in the color-singlet channel obtained with Eq. JSJ for 
the 6 different values of the temperature. 

singlet free energy 2(], 2?J can be reasonably well reproduced by the the following form reminiscent 



of a screened Cornell potential (as suggested, e.g., in Ref. 28]), 



F l (r,T) = -Ze-<*(r,T)r + - e ~^ T >) , (1) 

r (j,(r, 1 ) 



with a "screening mass" 



a 



Mr , r) = ^ e -0.3/r (2) 
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and two fitting functions given by 



a(r,T) 



b(t) = 1.1 - 3.6i 



4.3 1 2 + 17.5 1 3 



(3) 



where t = T/T c , a = 0.4 and a = 1.2 GeV 2 . The left panel of Fig. n summar i zes our fit to the 
lattice "data". Also shown is the unquenched zero-temperature potential as obtained in Ref. 
(recall that for T=0, E\ = F\, see also below), which is used to normalize the finite- T results 
at short distances, r < 0.2 fm, where the free energy is not expected to depend on temperature 
anymore. Our parametrization, Eq. (JJ), accommodates this T=0 constraint. 

As mentioned in the Introduction, the appropriate quantity in relation to the free energy that 
can serve as an effective potential appears to be the (color-singlet) internal energy E\. Following 
Kaczmarek et. al. we subtract the entropy contribution to the free energy according to 



F, - T- 



(4) 



The nonzero asymptotic value of the internal energy can now be interpreted as an in-medium quark 
mass that should not be included in the interaction part of the potential. One therefore assumes 
that the potential in the color-singlet channel can be extracted via 



V 1 {r,T)=E 1 {r,T)-E 1 (^,T). 



(5) 



The singlet potential is shown in the right panel of Fig. Q for the same values of temperature as 
the singlet free-energy (left panel). The potentials are appreciably larger in magnitude than the 
corresponding free energies and decrease with increasing temperature. 

To illustrate uncertainties in the determination of the potentials we compare in Fig. [5] our 
results with the ones obtained by other groups for temperatures of 1.5 T c (left panel) and 2 T c 
(right panel). While the potentials of Refs. 0] (Wo) and 0] (MP) are extracted from quenched 



1QCD, the one of Ref. [2JJ (SZ) and ours (MR) result from unquenched simulations. At 1.5 T c 
our potential is about 30-40% more attractive than SZ at distances between 0.1-0.8 fm, while 





FIG. 2: Comparison of Q-Q color-singlet potentials as determined from various lattice data in our (MR) 
and other works (SZ HJ, MP [3 and Wo [13); the left (right) panel is for T = 1.5 T c (2T C ). 
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similar deviations also occur between the quenched-based results, while at 2 T c our potential is 
less attractive than both the quenched (Wo) and unquenched (SZ) results. We therefore conclude 
that the current uncertainty in the extraction of the potentials amounts to about 50%, and that 
it does not yet allow for a systematic discrimination between quenched and unquenched results 
(provided the temperature dependence is normalized to T c , which is, of course, quite different in 
quenched (~260 MeV) and unquenched (~170 MeV) simulations). The pertinent uncertainties will 
be assessed below by performing T-matrix calculations for charmonium (c-c systems) with both 
our and the quenched Wo potential (cf. Sect. IIV Aj) . and comparing them to spectral functions 
from 1QCD jfj] obtained with different methods; reasonable agreement will be found. 

Toward a more complete description of the q-q interactions in the QGP we will in this work 
also consider the (repulsive) contributions from the color-octet channel. However, as pointed 
out in Ref. [3l]], the octet potential cannot be straightforwardly inferred from the Polyakov loop 
correlators. Due to a lack of better knowledge of the octet free energies, we here assume that the 
octet potential follows the leading-order result of perturbation theory, 

F 8 = -±Fi . (6) 

Again, we will check the sensitivity of our calculations to this approximation, by varying the 
coefficient in Eq. © by a factor of 0.5-2. 

For non-static quarks it is also important to include relativistic corrections |2fll . • Following 
Ref. 20 1 ] we implement a velocity- velocity interaction term by the replacement V(r) — * V(r)(l — 
&i ■ a.2) where a\ and 0.2 are quasiparticle velocity operators. As pointed out in Ref. [2l|, this 
procedure is strictly speaking correct only for a Coulomb-type potential. 



III. REDUCED BETHE-SALPETER EQUATION, QUARK SELF-ENERGY AND SELF- 
CONSISTENCY 

To evaluate quark-antiquark interactions in the QGP we employ the T-matrix approach, as is 
well known from the nuclear many-body problem. In relativistic field theory, the starting point is 
a 4-dimensional Bethe-Salpeter (BS) equation, 

T = K + J KSST , (7) 

where K denotes the interaction kernel and S is the single-particle propagator. Both quantities 
carry, in principle, dependencies on temperature and (baryon-) density of the surrounding medium. 
Since the effective q-q potential constructed in the previous section is essentially non-relativistic 
in nature, it is appropriate to employ the ladder approximation to Eq. (JJJ) in connection with 
neglecting virtual particle-antiparticle loops. Thus, we will identify the kernel K with the potential 
V with appropriate approximations in the propagator and scattering equation to be discussed in 
the following. 

The medium effects in the quark propagator, S, are encoded in a self-energy which we decompose 
according to 

^ = t + JtS. (8) 

The first term, S, represents a "gluon-induced" contribution due to interactions of (anti-) quarks 
with surrounding thermal gluons. In this work we do not calculate this term explicitly, but we 
will study how different (purely real) values affect our results. We note that a perturbative (hard- 
thermal-loop) form of this (mass-) term is widely used as a parameter in quasiparticle descriptions 
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FIG. 3: Schematic representation of the self-consistency problem composed of the Bethe-Salpeter equa- 
tion Q in ladder approximation and the quark self-energy, Eq. (J8J. Thick lines represent full fermionic 
propagators. The four different blocks correspond to the T-matrix (T), potential (V), self-energy (£) and 
"gluon- induced" self-energy (£). 



of the QGP EoS 0, 0, \M . The second term on the right side of Eq. (JSj) is the contribution 
to the self-energy induced by interactions with antiquarks of the heat bath which we compute at 
the same level of approximation as the T-matrix. In principle, the quark self-energy also receives 
contributions from interactions with thermal quarks (which could be significant especially in the 
scalar diquark channel), but we neglect them in this work. We also constrain ourselves to the 
case of vanishing quark chemical potential, fj, q = 0, which implies equal self-energies for quarks 
and antiquarks. The two equations ((7J) and (JHJ) constitute a self-consistency problem which is 
diagrammatically illustrated in Fig. |31 

Let us discuss the single-particle quantities in more detail. The full quark propagator obeys a 
Schwinger-Dyson equation, 

S = S + S ZS, (9) 

where So(k) = [|-rao] _1 . In the following, we set the current quark mass, mo, to zero and assume 
the self-energy to take the (chirally invariant) form 

T,(u, k) = a(u>, /c)7o + b{uj, k)h :• 7 , (10) 

since scalar and tensor contributions are suppressed due to chiral symmetry restoration (above T c ), 
whereas pseudoscalar and axialvector terms are absent due to parity invariance. The self-energy 
can be further decomposed as (cf., e.g., Ref. p^ ) 

7o S(^,k) = S + (^,A:)A + (k) - S_(^,A:)A_(k) , (11) 

where 

A ±( £) = i±^l (12, 



6 



are projectors on quark states with chirality equal (A + ) or opposite (A_) to their helicity, and 
Y>± = b(u>, k) ± a(uj, k). The quark propagator then follows as 



S(u, % = A+(cj, fe)A+(k) + A_(w, fc)A_(k) (13) 

with A-j- = — (a; =F (A; + E-i-)). Since the potential extracted from 1QCD is independent of chirality, 
the self-energy satisfies £+ = — that is b(LJ,k) = 0. Recalling Eq. ()10|) . this implies that 
the non-perturbative interactions only contribute to a chirally invariant (thermal) mass term for 
anti-/quarks. For the "gluon- induced" self-energy, £ in Eq. (j5J), we adopt a form suggested by 
the high-temperature hard-thermal loop result, characterized by a mass term, m, in the pertinent 
dispersion relation, 



\/k 2 + m 2 , (14) 



ignoring possible imaginary parts. Our default value for m is 0.1 GeV. 

Let us now turn to the scattering equation Q. As mentioned above, we neglect the (virtual) 
antiparticle components in the quark propagator and apply a 3-dimensional (3-D) reduction scheme 
to the 4-dimensional BS equation, facilitating its numerical evaluation substantially. The resulting 
Lippmann-Schwinger equation takes the form 

f d 3 k 

T a (E- q', q) = V a (q[ , q) - J V a (q! , k) G qq (E; k) T a (E; k, q) [1 - 2f(u k )} , (15) 

where E denotes the center-of-mass {CM) energy and q and q' are the in- and outgoing (off-shell) 
3-momenta in the CM (as usual, the on-shell T-matrix is defined by q=q' with E = 2uj q where uj q 
is the on-shell single-quark energy); a = 1,8 labels color-singlet and -octet channels, and 

/M = ^TI (16) 

is the Fermi-Dirac distribution. The explicit form of the two-particle propagator, G qq -(E;k), de- 
pends on the 3-D reduction scheme. Unless otherwise stated, we adopt the Blankenbecler-Sugar 
(BbS) [3] prescription leading to 

G qq -(E;k) = Uk — (BbS) , (17) 

but we have checked that our results are very similar when employing the Thompson scheme |3_ 
with 

G qq (E;k) = l *. v -, 7T (Th). (18) 

2 oj k - E/2 + zS/(wfc, k) 

In both Eqs. (|17|) and ()18j) tu k denotes the on-shell quasiparticle dispersion law, i.e., the solution 
of the equation 



Vk 2 + m 2 + Z R (oj k ,k), (19) 



with Y,r and T,j the real and imaginary part of the self-energy. Finally, the potential figuring into 
Eq. ()15|) follows from our 1QCD parametrization via Fourier transformation, 

y a (q',q)= / d 3 rV a (r)e^-*'> r . (20) 
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To solve Eq. (|15|) it is convenient to work in a partial-wave basis. Expanding T-matrix and 
potential, 

V a (q!,q) = 4vr£(2Z + 1) V a j(q',q) P,(q' • q) , (21) 
I 

r«(£;q',q) = 4tt]T(2Z + 1) T ajl (E;q',q) P(q'-q) , (22) 
I 

allows to perform the angular integrations to yield 

T a i{E; q, q) = V a ,i{q , q) - - [ k 2 dk V ajl (q',k) G q q{E;k) T a>l (E;k,q) [1 - 2f(u k )] . (23) 

vr J 

In the present study we will constrain ourselves to 5-wave channels, deferring higher waves to 
future work. In Ref. [21] it was found that P-wave bound state formation is strongly suppressed 
(in accordance with our own estimates). Concerning spin-isospin channels, we recall that in the 
chirally restored phase the spectral functions of chiral partners {e.g ., n-a, p-a\) degenerate, which 
is also reflected in the spectral functions extracted from 1QCD |4Q|. Within the naive constituent 
quark model, ir and p states are S-w&ve q-q bound states, whereas a and a\ are in a P-wave state. 
Interestingly, 1QCD spectral functions find an additional (approximate) degeneration of ir and p 
states above T c In as far as an interpretation of these objects as qq states applies, this might 

be taken as an indication for a spin-symmetry much like in heavy-quark effective theories. In view 
of these considerations, and due to the fact that our lQCD-extracted potential is flavor-blind, we 
will assume the color-singlet S-wave states to appear with a spin-isospin degeneracy corresponding 
to TT+p states, dsj = 12. Since the color-octet potential does not carry any flavor-dependence 
either, the same factor will be applied to the color-octet states. 

With the q-q T-matrix at hand, we can proceed to calculate the explicit expression for the quark 
self-energy due to interactions with anti-quarks. Within the imaginary time formalism the latter 
follows from closing the forward scattering T-matrix with a thermal q propagator, 

= jfdaj ^ (-T) £T^ + ^;p,p') D,(v,p') (24) 

(here, z u = iri(2v + 1)T are fermionic Matsubara frequencies, d\% = 1,8 is the color degeneracy 
factor and the factor 1/12 represents the average over the 3x2x2 (color x flavor x spin) initial 
quark states). Using the spectral representations of both T-matrix and q propagator to perform 
the Matsubara sum, and after analytic continuation to the real axis, the self-energy takes the form 

E a ( W; p) = ^ da r ±l r d A r A{u/t fc) m±m lmTaiE . k + p) (25) 

V F> 12 J 2vr J vr J (2tt) 3 y J to + to' - E + irj V V ' V ' 



with the Bose distribution 



and the quark spectral function 



9W = jW^i (26) 



A( U ,k) = -2^,k) 

V ' (u - Vk 2 + m 2 - V R (u, k)) 2 + Ej(w, A;) 2 V ' 

To further simplify our task we assume in the following a quasiparticle approximation for the 
spectral function, 

A(u,k) = 2tt5(uj -oo k ), (28) 
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where Uk is obtained from the self-consistent solution of Eq. ()19|) (we will check this approximation 
below). If we furthermore neglect the (weak) energy dependence of g(E) close to the pole of the 
principal value integral in Eq. (|25|). we can recover the real part of T q q to cast the self-energy in 
compact form, 

X a (u;p) = j±d a J Lf(o, fc ) + g( u + T a qq (E) , (29) 

where x = cos 9 (with = Z(p, k)) and the CM energy of the on-shell T-matrix is given by 

E = yV fc + w ) 2 - (p + k)2. (30) 



IV. T-MATRIX, SELF-ENERGY AND SPECTRAL FUNCTION 



In this section we discuss the numerical solutions to the set of equations l|19|). (|23jl and (|29|). 
Self-consistency is achieved by iteration, starting with the calculation of the T-matrix using a 
constant self-energy in the first step. The self-energy is then calculated from ((2*9"|) and used to solve 
the on-shell condition (|19|) . The pertinent quasiparticle dispersion-law is then re-inserted into the 
T-matrix equation and the procedure is iterated until T-matrix and self-energy converge (typically 
within less than 10 iteration steps; we have also verified that the final results are insensitive to the 
initial input value for the self-energy). 



A. Quark- Antiquark T-matrix 

The T-matrix equation (|23j) is solved using the matrix inversion algorithm of Haftel and 
Tabakin [4]J (after discretizing the momentum integration). To assess the possible forma- 
tion of bound states, the T-matrix needs to be calculated below the nominal q-q threshold, 
E t hr = 2(m + T,ji(Ethr/^, 0)). The potential does not depend on the CM energy E, and, due 
to its nonrelativistic character, is only defined for real external 3-momenta q and q' ■ We therefore 
define the subthreshold on-shell T-matrix by setting the external momenta q=q'=0. In the follow- 
ing we will refer to a peak in the imaginary part of the T-matrix as a bound-state (resonance) if 
the energy of the maximum is located below (above) the quasiparticle threshold, Ethr- 



1. Charmonium Systems 



To check the reliability of the parametrization of the potential in the singlet channel, and of the 
algorithm to compute the T-matrix, we first apply our approach to the c-c (charmonium) sector 
by using a (constant) quark mass of m = 1.8 GeV which approximately reproduces the vacuum 
J/iJj mass at the lowest temperature (note that self-consistency does not play a role here since 
the thermal abundance of c-quarks is strongly suppressed; for numerical purposes, we used a fixed 
imaginary value for the self-energy, £/ = —10 MeV, and = 0). The results are displayed in Fig.|I] 
for three different temperatures, 1.2 T c (left panel), 1.5 T c (middle panel) and 2T C (right panel). 
As the temperature increases the charmonium state moves up in energy (reflecting a decreasing 
binding energy) reaching the threshold (Ethr = 3.6 GeV) at T ~ 2T C after which the resonance 
peak essentially dissolves (also note that the strength in the T-matrix is much reduced at 2 T c as 
compared to the lower temperatures). This behavior is in reasonable (qualitative) agreement with 
both 1QCD calculations 0] and effective potential models using a Schrodinger equation 17. 
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2.5 3 3.5 4 2.5 3 3.5 4 2.5 3 3.5 4 

E [GeV] E [GeV] E [GeV] 



FIG. 4: Real (full red line) and imaginary part (absolute value, dashed blue line) of T-matrix in the color- 
singlet channel for charmonium (with a charm-quark mass of m — 1.8 GeV) at T = 1.2 T c , T = 1.5 T c , 
T = 2T C (left, middle and right panel, respectively) as a function of CM energy E, based on our potential 
paramctrization extracted from unquenched 1QCD. 




2.5 3 3.5 4 2.5 3 3.5 4 ~ zu 2.5 3 3.5 4 

E [GeV] E [GeV] E [GeV] 



FIG. 5: Real (full red line) and imaginary part (absolute value, dashed blue line) of T-matrix in the color- 
singlet channel for charmonium (with a charm-quark mass of m = 1.8 GeV) at T = 1.2 T c , T — 1.5 TL, 
T = 2 T c (left, middle and right panel, respectively) as a function of CM energy E, based on a potential [lj| 
extracted from quenched 1QCD. 



To check the sensitivity to the underlying potential (recall the discussion in Sec. [Q] around 
Fig. I2J), we have repeated the calculations for the charmonium T-matrix in the singlet channel using 
the (quenched-based) potential of Ref. 0] , cf. Fig.El At the lower temperature of 1.2 T c the binding 
is significantly less pronounced (by about 0.25 GeV) as compared to our parametrization, as to be 
expected from the less attractive potential. At higher temperatures the agreement improves, and 
both potentials lead to a very similar temperature where the state crosses the c-c threshold (close to 
2T C ), with strongly reduced strength. This, in turn, is again in line with the Schrodinger-equation 
approach, in which the c-c system becomes unbound around ~2T C 17]. While the resonance at 
2 T c appears to be rather narrow, we recall that we did not include here (temperature dependent) 
absorptive parts and reduced masses for the c-quarks (nor inelastic charmonium reaction 
channels |42j), all of which are expected to increase the width of the charmonium states. 
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0.5 1 1.5 2 " iu 0.5 1 1.5 2 0.5 1 1.5 



E [GeV] E [GeV] E [GeV] 

FIG. 6: Real (full red line) and (absolute value of the) imaginary part (dashed blue line) of the light-quark 
(on-shell) T-matrix in the color-singlet channel at temperatures T = 1.2 T c , T = 1.5 T c and T = 1.75 T c 
(left, middle and right panel, respectively) as a function of the qq CM energy E, with a "gluon-induced" 
quark-mass term to = 0.1 GeV. 
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FIG. 7: Light-quark T-matrix in the color-octet channel vs. qq CM energy at T — 1.2 T c , T — 1.5 T c 
and T = 1.75 T c (left, middle and right panel, respectively) with to = 0.1 GeV. Solid (red) line: real part; 
dashed (blue) line: imaginary part (absolute value). 



2. Light- Quark Systems 

Turning to the light-quark sector, the self-consistent results for real and imaginary part of 
the on-shell T-matrix for quasiparticles with a gluon-induced mass-term of m = 0.1 GeV are 
summarized in Figs. El and [7| for temperatures T = 1.2 T c , T = 1.5 T c and T = 1.75 T c . 

At T = 1.2 T c the color-singlet T-matrix exhibits a relatively narrow bound state located signif- 
icantly below the q-q threshold energy of E t h r ~ 0.52 GeV (corresponding to twice the real part of 
the total quark self-energy discussed below). When increasing the temperature to 1.5 T c , the state 
moves to higher CM energy above the threshold {Ethr — 0.48 GeV) which, not surprisingly, is 
accompanied by a significant broadening. Note also that the peak value is substantially reduced as 
compared to the 1.2 T c case, substantially more than to be expected from the broadening alone. We 
assign this behavior to the decrease in the potential, cf. right panel of Fig. ^ reflecting an overall 
reduction in interaction strength. The trends in suppression, broadening and upward energy-shift 
continue at T = 1.75 T c where the resonance has now essentially melted as indicated by a width of 
almost 1 GeV, comparable to its mass. These results may be put into context with computations of 
mesonic spectral functions in (quenched) lattice QCD. For (reasonably) light quarks jsEO]) their 
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main features above T c are a gradual increase of the peak position (corresponding to the "meson 
mass") with temperature (roughly proportional to T), accompanied by a broadening. The bound 
/ resonance states depicted in Fig. approximately share both of these features. 

The T-matrix in the color-octet channel is displayed in Fig. [7| for the same set of temperatures. 
As to be expected for a purely repulsive potential, we find a smooth (non-resonant) dependence of 
both real and imaginary part with CM energy (with a substantial suppression at higher T, as in 
the singlet case). The imaginary part is very small, and also the real part appears to be small when 
compared to the singlet channel. We recall, however, that the octet contribution to the self-energy, 
Eq. (|29j). enters with a weight which is by a factor of 8 larger than for the singlet one, rendering it 
an important effect as will be seen below. 



B. Self-Energy 



We proceed to the single-quark self-energies as calculated from the interactions with antiquarks 
of the heat bath using the expression, Eq. (|29[). based on the self-consistent 5-wave q-q T-matrices 
in the "7r" and "p" channels as obtained in the previous section. We recall that the real part of the 
self-energy corresponds to a chirally invariant mass-term, whereas its imaginary part determines 
the width of a quark (-quasiparticle) according to T = — 2 ImS. Since we work at zero quark- 
chemical potential, n q = 0, the same results hold for antiquarks. We also recall that our on-shell 
approximation scheme for the self-energy implies that the effects of bound states are not captured 
by Eq. (|29j) . since in the integration over the T-matrix only energies above the q-q threshold, Ethr, 
contribute. 
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FIG. 8: Real (solid line, red) and imaginary (dashed line, blue) part of the on-shell quark self-energy as a 
function of 3-momentum at temperatures T = 1.2 T c , T — 1.5 T c and T = 1.75 T c (left, middle and right 
panel, respectively) with m = 0.1 GeV. 



In Fig. |H] the on-shell self-energy is displayed for the same selection of temperatures as in the 
previous section. Both real and imaginary part are smooth functions of the quark 3-momentum with 
maximal values at k = 0. Note that real part is positive, implying that the repulsive contribution 
from the octet channels overcomes the attractive singlet channels. The imaginary part (width), on 
the other hand, chiefly arises due to resonant scattering in the singlet channel. 

More quantitatively, in the temperature regime 1.2-1.5 T c , the nonperturbative contribution 
to the thermal quark mass reaches values of around 150 MeV at small momenta, decreasing to 
~50 MeV at 1.75 T c . With the underlying "gluon-induced" mass term of m = 100 MeV, the total 
thermal mass, m + amounts to 150-250 MeV. This is smaller than effective (perturbative) 
thermal quark masses required in phenomenological fits to the QGP EoS of 1QCD jia. liil liil liil] . 
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To improve upon this, we have performed self-consistent calculations with a gluon- induced mass 
term of m = 250 MeV. It turns out that, at given temperature, the "mesonic" states are slightly 
stronger bound, but in general the behavior of T-matrix and pertinent self-energy are quite similar 
to those obtained with m = 100 MeV. E.g., for T = 1.5 T c (cf. Fig. 03J), the resonance structure is 
right at threshold, the quark width reaches almost 200 MeV, and the combined real part at low 
momenta amounts to a quark mass of m + ~ 350 MeV. 

An important aspect of our results are the rather large imaginary parts of the quark self-energy, 
translating into widths of about 200 MeV at low momenta for temperatures around 1.5 T c . As 
mentioned above, the width is almost entirely generated by the resonant scattering in the singlet 
channel; this is nicely illustrated by the significant increase in ImS when going from 1.2 to 1.5 T c (cf. 
left and middle panel in Fig.EJ), during which the state in the T-matrix moves from below to above 
threshold (cf. left and middle panel in Fig. |SJ), i.e., converts from bound state to resonance 1 . The 
magnitude of the quark widths is quite comparable to the thermal masses, qualitatively supporting 
the notion that the QGP could be in a liquid-like regime [43I. Even at the highest considered 
temperature of 1.75 T c , and at typical thermal momenta (k ~ 3T ~ 0.9 GeV), the quark width 
due to scattering off antiquarks is between 50 and 100 MeV. This does neither include P-wave 
interactions, nor strange antiquarks, nor any contributions from scattering off quarks or gluons. 

Finally let us come back to the uncertainty associated with the interaction in the octet channel 
related to the perturbative ansatz, Eq. @. If the coeffiecient in Eq. © is increased (decreased) 
by a factor 2 (for T = 1.5 T c and m = 0.25 GeV), the imaginary part of the self-energy barely 
changes, whereas its real part increases (decreases) by about a 40%. 
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FIG. 9: Real (full line, red) and imaginary (dashed line, blue) parts of the T-matrix in the color-singlet 
channel (left panel), color-octet channel (central panel) and corresponding (singlet+octet) self-energy (right 
panel) at a temperature T = 1.5 T c using a "gluon-induced" mass term of m — 0.25 GeV. 



C. Quark Spectral Functions and Normalization Condition 

To better elucidate the validity of the quasiparticle approximation, Eq. (|28|). we compute the 
off-shell real and imaginary parts of the self-energy using Eq. (|29|) and obtain the pertinent quark 
spectral function A(co, k) from Eq. (|27|) . In Fig. ^3we depict A(tu, k) as a function of quark energy, 



1 We recall that bound states are not accessible in on-shell 2^2 scattering; even if a resonance is close to threshold 
it does not contribute effectively to rescattering processes if the average thermal energy of particles from the heat 
bath is significant. The contribution of bound states to the self-energy can be included rather by going beyond 
the quasiparticle approximation, i.e., evaluating Eq. 12511 with the off-shell spectral function, Eq. 1271 . 
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FIG. 10: Off-shell spectral function A(ui, k) as given by Eq. (|27|) vs. quark energy for different values of the 
quark momentum at T = 1.5 T c and for m = 0.25 GeV. 

u, for various fixed momenta at a temperature of T = 1.5 T c . On the one hand, this reiterates 
the large effect of the width for low momenta and calls for an off-shell treatment to improve the 
reliability of our results in the (sub-) threshold region of the T-matrix. On the other hand, for 
larger momenta (including typical thermal momenta) the quasiparticle approximation as applied 
in our calculations appears to be reasonably well justified. 

As another check of our approximations we have evaluated the norm of the quark spectral 
functions defined by 



The unitarity condition for A(u, k) requires I(k) = 1 for each momentum k. This relation is rather 
well satisfied, I(k) > 94 %, for all momenta considered in Fig. 

V. CONCLUSIONS AND OUTLOOK 

In the present article we have set up a self-consistent many-body scheme of Brueckner-type 
to assess nonpertubative properties of (anti-) quarks and mesonic composites in a Quark-Gluon 
Plasma at temperatures T ~ 1.2-2 T c . Our key ingredient to describe the q-q interaction in the 
QGP was a driving kernel (potential) extracted from unquenched finite-T lattice QCD calculations 
for the free energy of a heavy-quark pair, supplemented with corrections for relativistic motion. 
Our main objective was to go beyond earlier applications to bound states by solving the scattering 
problem thereby accounting for absorptive effects (finite imaginary parts). The self-consistent set 
of single-quark Dyson and two-body scattering equations has been solved by numerical iteration 
employing a nonrelativistic reduction of the Bethe-Salpeter equation in connection with a quasi- 
particle approximation for the quark propagators. One of our main new findings is that the 1QCD 
potentials (dynamically) generate S'-wave resonance states above the q-q threshold up to tempera- 
tures of ~ 2 T c . These resonances (assumed to occur with a degeneracy corresponding to 'V- and 
"//'-mesons), in turn, play a key role in inducing large quark scattering rates (=imaginary parts 




(31) 
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of the quark self-energy) as indicated by single-particle widths of V ~ 200 MeV at temperatures 
around 1.5 T c . At the same time, significant (positive) real parts arise from repulsive interactions 
in the color-octet channel entailing thermal masses of up to ~150 MeV. We expect that additional 
contributions to the quark mass of ~250 MeV (induced, e.g., by interactions with thermal gluons 
as parametrized in quasiparticle models) will be necessary to account for the QGP EoS computed 
in lattice QCD. Nevertheless, especially at low momenta, the quark widths are comparable to 
the thermal masses, which could be indicative for liquid-like properties of the QGP at moderate 
temperatures. 

Our analysis suggests several directions for future work. First, the accuracy of our approxima- 
tions should be scrutinized. This includes improving upon the quasiparticle approximation of the 
quark spectral function by implementing its off-shell (energy-) dependence (as, e.g., carried out in 
Ref. [45j for a hot pion gas), most notably at low energies to incorporate bound-state contributions 
to the quark self-energy. The scattering equation ought to be extended to finite total 3-momentum 
of the mesonic composites. Even though we expect the q-q channel to constitute a major part of 
the in-medium interaction, a more complete treatment including q-q and q-g channels is desirable. 
It is also conceivable that processes of the type qq — > Mg (inverse gluon-dissociation; M: mesonic 
state) could be significant, as they render bound states accessible in (on-shell) 2-body scattering. 
In a broader context, the underlying EoS of the interacting system needs to be investigated, which 
is obviously not an easy task. On the phenomenological side, to address the problem of early 
equilibration at RHIC, it will be of great interest to calculate the thermal equilibration timescales 
for (anti-) quarks based on the resonant scattering amplitudes found here (e.g., within a Fokker- 
Planck equation). The elastic scattering rates of around l/(fm)/c as found in this work, together 
with the isotropic angular dependence inherent in 5"-wave rescattering, look promising. For gluons 
the situation could be more involved since, besides bound states as suggested in Ref. other 
thermalization mechanisms might be operative, e.g., gg <-> ggg processes 4f], 4^|. In this respect, 



charm quarks are of particular importance, as their number is presumably frozen after primordial 
production, and genuine 2— >3 processes are absent. Indeed, the recent analysis of Ref. [lfj] has 
shown that "D" -meson resonances in the QGP can accelerate thermal relaxation times obtained 
from pQCD by a factor of ~3. A rather straightforward extension of our approach to the heavy- 
light sector should therefore be pursued. The formation of mesonic composites in the cooling QGP 
phase of a heavy-ion collision could furthermore serve as a "pre-hadronization" mechanism, and 
thus improve phenomenologically successful quark-coalescence models at RHIC 0, 0, (e.g., 
with respect to the question of energy conservation). Significant future efforts will be required 
to possibly develop such a scheme into a quantitative phenomenology. Further progress will also 
reside on increasing information from finite-T lattice QCD to provide both input and constraints 
to a many-body approach as presented here. Clearly, a thorough understanding of the intricate 
properties of the strongly interacting matter above T c , and its implications for ultrarelativistic 
heavy-ion experiments, is an exciting future task. 
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